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Abstract. The stages that follow the merging of two neu- 
tron stars are discussed. It is shown that if a rapidly ro- 
tating gravitationally bound object is formed after the 
merging (a spinar or a massive neutron star), then the 
characteristic time of its evolution is determined by a fun- 
damental value 



7-10^ 



where the dimensionless value n depends on the exact 
equation of state of nuclear matter. The hypothesis is dis- 
cussed as to whether the residual optical emission of the 
gamma-ray bursts is pulsar-like and its evolution driven 
by magnetodipole energy losses. It is shown that binary 
neutron star mergings can be accompanied by two gravita- 
tional wave burst separated either by the time of spinar's 
collapse tspin or neutron star cooling time (~ 10 s), de- 
pending on the masses of neutron stars. 
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1. Introduction 

The detection of optical and X-ray emission after the 
gamma ray bursts GRB 970228, GRB 970508 (Groot et 
al. 1997a; Groot et al. 1997b; Metzger et al. 1997b ; Costa 
et al. 1997b; Sahu et al. 1997; Bond 1997; Galama et 
al. 1997; Djorgovski et al. 1997a; Metzger et al. 1997c; 
Schaefer et al. 1997; Djorgovski et al. 1997b; Djorgovski 
et al. 1997c; Groot et al. 1997; Donahue et al., 1997) may 
be interpreted in terms of the formation of a transient 
rapidly rotating gravitationally bound object — a heavy 



neutron star (NS) or spinar — an object with the equilib- 
rium maintained either by the fast rotation ("cool" spinar, 
CSP) or by both rotation and pressure ("hot" spinar, 
HSP). 

Let us assume that two neutron stars with masses Mi 
and M2 are merging. The following state of the after- 
merging object is determined by the ratio of the result- 
ing total mass and the Oppenheimer-Volkoff limit. Two 
different scenarios may be envisaged as follows: 



Ml + M2 > Mov 
Ml + M2 < Mov 



(A) 
(B) 



Here and below we interpret the Oppenheimer-Volkoff 
limit not as the standard value derived for the cold equa- 
tion of state of baryonic matter for a non-rotating neu- 
tron star, but as a modified one. In the general case the 
Oppenheimer-Volkoff limit is a function of the angular 
spin velocity of the object, its entropy, and the specific 
equation of state: Mqv = Movi^^iS^EqSt). 

Each neutron star can have a mass lying between the 
limits: 



M„ 



< Mi,M2 < Mov 



Send offprint requests to: G.V. Lipunova 



The value of Mmm ~ O.2Af0 was derived by Landau 
(1938). In a standard modern scenario, it is commonly sug- 
gested that Mmin ~ 1.2M0. Thus we can expect the differ- 
ent evolutionary tracks depending on the specific masses 
of NS. 

2. Mergingology 

2.1. Case (A) 

In this case, we can expect that after the merging a black 
hole results from a direct collapse during the time ~ 10~^s 
and that the most energy is emitted in the gravitational 
wave burst. This scenario is discussed more frequently in 
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the literature, and GRB phenomenon can be related with 
the relativistic particle ejection in the form of a Fireball 
(Rees & Meszaros, 1992) or a beam of protons (Shaviv & 
Dar, 1996). In addition, a certain fraction of radiated en- 
ergy can be related with the pulsar mechanism (Lipunov 
& Panchenko 1996; Lipunova 1997). No gravitationally 
bound object can be formed in this case outside the hori- 
zon. We can present these stages by the following way: 

NS + NS ^ BH + GWB + GRB + vB 

{GWB - gravitational waves burst; vB - neutrino burst). 

From our point of view the more interesting scenario 
is the Case (B). 

2.2. Case (B): Mi + M2 < Mqv 

This variant can be realized if cither two merging neutron 
stars have small masses or the Oppenheimer-Volkoff limit 
is very large. 

Is it possible that the Oppenhcimcr -Volkoff limit ex- 
ceeds 3M0? First, it is known (Friedman & Ipser, 1987) 
that the fast rotation (which is naturally expected after 
the merging) increases the Oppenhcimcr Volkoff limit to 
the value ~ SM© for hard equations of state. Second, the 
object formed is not degenerate due to its high tempera- 
ture and the equilibrium is maintained both by fast rota- 
tion and entropy ("hot" spinar). And last, Oppenheimer- 
Volkoff limit can be high because of relativistic behavior 
of nuclear forces. 

Thus we can present these three sub-scenarios as fol- 
lows: 



Finally, there is a case of a high Oppenheimer-Volkoff 
limit for the cool non-rotating object. 



Mi+M'2< Mov 



always ! 



This variant leads to the formation of a very powerful 
pulsar (maybe without pulsation) with the maximum spin 
rotation. 



NS + NS- 
The characteristic time 



PSR 



i-spm of its evolution is gov- 
erned by the momentum loss rate. 



3. The rate of the angular momentum losses 

In both cases of a cool spinar (Lipunova 1997) and of a 
fast-rotating NS, the specific time of their evolution is 
determined by the rate of magnetodipole energy loss 



diu] 



and 



ts 



2dj 



■3 C3 



2 Mc^ 



5 BlRioj'^ ■ 



We assume: 

the inertia moment 
the mass 



I = f Mi?2 , 
M = Mi+ M2, 
the magneto-dipole moment fi = BoR^/2 . 
The angular spin velocity of the post-merging object must 
be close to the limit: 



NS + NS^ HSP + GWB + GRB + vB 

NS + NS^ CSP + GWB + GRB + vB 

NS + NS ^ NS + GWB + GRB + vB 

Let us consider the case of the HSP. Its lifetime is com- 
pletely determined by the cooling time which, according to 
different calculations, is of the order of ~ 10 s (Shapiro & 
Teukolsky, 1983). Then, in the time interval tcooh the col- 
lapse accompanied by the GWB, neutrino emission, and 
possible weak photon emission can be expected: 

HSP BH + GWB + vB + -i 

It seems very attractive to identify this cooling time 
with the mean characteristic gamma-ray burst duration 
~ 1 ^ 10 s! 

Second, the most interesting sub-scenario is when the 
centrifugal forces make the main contribution to the equi- 
librium ("cool" spinar). In this case the lifetime of the 
spinar is completely defined by the characteristic time of 
the angular momentum loss tspin and evolutionary track 
looks like 

CSP -^BH + GWB + -f + e+ + e- +iy 



CO = {GM/Rl) 



.3^1/2 



Then we obtain: 



t. 
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spin 



5BiGR 



2 • 10^ 



B 



\ lO^CTO 



Thus, this duration is determined mainly by the mag- 
netic field. If we assume that a gravitationally bound ob- 
ject magnetic field is equal to the critical value close to 
the Schwinger limit: 



Ti 



eB, 



rrifC 



Bf. 



4.3 • 10^^ G . 



Expressing the radius and the mass of the NS in terms 
of fundamental constants we obtain the fundamental value 
for the lifetime of such an object: 



mtd'/^Gy^ 



7.6 • lO^s . 



Taking into account the real mass of NSs and specific 
equation of state, this time can be modified as 



ts 



T-K , 
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Spinar 

NS spin down 

1.2 I 1 1 1 1 1 1 1 r 




t, days 

Fig. 1. The spinar luminosity evolution. The magnetic field 
B = 4.3 • 10" Gs at i? = 20 km. 



where k depends on the exact equation of state of nuclear 
matter. This duration accords with the specific fundamen- 
tal value of luminosity. 

4. GRB light curve 

Here we present the alternative model to the now fre- 
quently discussed model of GRB — the radiation of a 
Fireball (see Meszaros & Rees 1997). We admit that the 
models of fast-rotating pulsar or spinar with extremely 
high magnetic fields do not wholly substitute the model 
of Fireball (especially, concerning the gamma burst itself) 
but accompany the process of radiation and, possibly, at 
later stages of GRB afterglow, dominate in a GRB spec- 
trum. Note, that these mechanisms can supply the emis- 
sion in a wide range of wavelengths, as radiopulsar studies 
confirm. 

We suggest that part of the observed optical and X-ray 
afterglow of a GRB can relate to the pulsar mechanism. 
As an emission from Fireball decreases to the undetectable 
level, the pulsar mechanism can become the main contri- 
bution to the afterglow. 

We can construct the luminosity evolution for a cool 
spinar collapse (Lipunova, 1997) and for a pulsar spin 
down. 

Supposing that the optical emission is produced by the 
pulsar mechanism acting with the critical magnetic field, 
one can derive: 

L^2- lO^^erg/s ( ^) RlP^j K{t, u) x 



I ' 1 

-10 10 20 30 40 50 

t, days 

Fig. 2. The neutron star spindown rate and luminosity evolu- 
tion. R = 30 km, Umax = 660 Hz. 



where Rq = {R/lO^cm), P1.5 = (P/1.5tos), 
Af3 = (M/3M0). The coefficient K{t,iy) is the ratio of 
optical radiation to the total energy loss by a pulsar. Of 
course, it is hard to expect the ratio of optical radiation 
to the total rotational energy loss to be constant, as evi- 
denced by radiopulsar studies. As it is, the real power of 
time dependence can vary from —2. 

Fig. 1 shows the characteristic times of luminosity de- 
creasing to be in a rather good correlation with the ob- 
served ones (see Groot et al. 1997a; Groot et al. 1997b; 
Metzger et al. 1997b Costa et al. 1997b; Sahu et al. 1997; 
Bond 1997; Galama et al. 1997; Djorgovski et al. 1997a; 
Metzger et al. 1997c; Schaefer et al. 1997; Djorgovski et 
al. 1997b; Djorgovski et al. 1997c; Groot et al. 1997; Don- 
ahue et al., 1997). The model of a neutron star spin down 
is calculated for the initial angular velocity ly — 660 Hz, 
which corresponds to the minimum spin period observed 
in millisecond pulsars. 

The lack of optical counterparts to other GRBs may 
be explained by another relation between the total mass of 
the system before merging and the Oppenheimer-Volkoff 
limit and, as a result, by another scenario of neutron star 
coalescence. 
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